Charged Multiplicities at SPS and RHIC and consequences for J/ijj suppression 
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Hadron multiplicities in nucleus-nucleus interactions are calculated in the Dual Parton Model 
and its dependence on the number of collisions and the number of participants is analyzed. 
Shadowing corrections are calculated as a function of impact parameter and the multiplicity 
per participant as a funtion of centrality is found to be in agreement with experiment at 
SPS and RHIC energies. The obtained results are used to compute the J/tp suppression in a 
comover approach. 



1 Charged Multiplcities in the Dual Parton Model 

In the Dual Parton Model (DPM) the charged multiplicity per unit rapidity in a symmetric 



Here P and T stand for the projectile and target nuclei, n(b) is the average number of binary 
collisions and the average number of participants of nucleus A. These quantities can 

be computed in a Glauber model, k is the average number of inelastic collisions in pp and 
//(6) = kn{b) /nA{h) is the total average number of collisions suffered by each nucleon. 

In the DPM each inelastic collision leads to two strings and the total number of strings 
is 2kn. The charged multiplicity produced by a single string is obtained by a convolution of 
momentum distribution functions and fragmentation function, and it depends on /x due to energy 
conservation. We see from (1) that the multiplicity is given by a linear combination of and 
n with coeficients which depend on the impact parameter b (via n{b)). 
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Figure 1: a): The values of dN'"^/drjJ ripart versus Upart for PbPb collisions at ^/s — 17.3 GeV in the range 
—0.5 < 77cm < 0.5 computed horn eq. (0), compared with the WA98 datau. b): The values of dN'^'' /drj / (0.5npart) 
for Au-Au collisions at ^/s = 130 GeV in the range —0.35 < r^cm < 0.35 computed from eqiri (hi) taking into 
account shadowing corrections (see main text). They are compared to the datau. 



Shadowing corrections in Gribov theory are universal!!, i.e. they apply both to soft and hard 
processes. They are closely related to the size of diffractive production and, thus, are controlled 
by triple Pomeron diagrams. The reduction of the multiplicity resulting from shadowing correc- 
tions has been computed ini. These corrections are negligeable at SPS energies but at RHIC 
energies they reduce the multiplicity by 40 to 50%. 

We present the resultsll obtained at two different energies : ^/s = 17.3 and 130 GeV. The 
corresponding non-diffractive cross-sections are a^o = 26 and 33 mb, respectively. We take 
k = 1.4 and 2.0 corresponding to dN^^^ /dy = 1.56 and 2.72. The result in absence of shadowing 
at ^/s = 17.3 is shown in Fig 1 (a) (solid line). We obtain a mild increase of the multiplicity per 
participant consistent with the results of the WA98 Collaboration!. This increase gets stronger 
with increasing energies. As we pointed out before, shadowing corrections are negligeable at SPS 
energies but their effect is large at RHIC. Unfortunately, shadowing corrections have a rather 
large uncertainty at RHIC energies. Two alternative calculations a of shadowing lead to the 
results at ^/s = 130 GeV shown by the solid lines in Fig 1 (b). Clearly, with the larger values 
of the shadowing corrections we obtain a quantitative agreement with the PHENIX dataB. 

Note that our calculations refer to dN/dy while the first RHIC measurements refer to dN/drj. 
The latter is, of course, smaller at mid rapidities. This difference is negligibly small as SPS where 
the laboratory pseudo-rapidity variable is used. However, at ^/s = 130 GeV where rjcm is used 
instead, their ratio can be as large as 1.311. 



2 J/ip suppression 

Now we are going to use these multiplicities, computed in the framework of DPM, to determine 
the J /ip suppression in the comovers approach. In this model, the J /ip survival probability is 
the product of two factors Sabs{b,s) ■ Sco{b,s). The first factor represents the suppression due 
to nuclear absorption of the cc pair. Its expression, given by the probabilistic Glauber model, 
is well known. It contains a parameter, the absorptive cross-section (Jabs- The second factor 
jS6p(6, s) represents the suppression resulting from the interaction with comovers. Its expression 
depends on the averaged interaction cross-section aco and on the value of the multiplicity at 



each b in the region of the dimuon trigger (0 < y* < 1), for which we take the results of section 
1, {Ny° is the charged multiphcity; the factor 3/2 takes care of the neutrals) 
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Experimentally, the ratio of J/ip over DY is plotted as a function of either Ex or the energy of 
the zero degree calorimeter Ezdc- Et is the transverse energy of neutrals deposited in the NA50 
calorimeter, located in the backward hemisphere (1.1 < yiab < 2.3). Using the proportionality 
between Et and multiplicity, we have 
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Here the multiplicity of comovers is determined in the rapidity region of the NA50 calorimeter. 
The energy of the zero degree calorimeter is defined as 



EzDcib) = [A- nA{b)]Ein + anA{b)Ein , (4) 

Here is the number of participants, A — ua the number of spectators, and Ein = 158 GeV 
is the beam energy. The last term represents the small fraction of wounded nucleons and/or fast 
secondaries that hit the ZD Calorimeter. The values of q and a are obtained, respectively, from 
a fit of the tails of the Et and Ezdc distributions. With the obtained values, the measured 
Et — Ezdc correlation is well reproduced. 

We see from Eqs. (^) and that, in order to describe the centrality dependence of the 
J /'4> suppression, it is paramount to have a good description of the b dependence of Ny" - both 
in the rapidity region of the dimuon trigger and in the one of the Et calorimeter. In the latter 
we obtain an scaling in the number of participants, i.e. we recover the wounded nucleon model. 

The model allows to compute the ratio J/ ip over DY versus either Et or Ez dc from periph- 
eral collisions up to the knee of the Et or Ezdc distributions. To go beyond it, we have to intro- 
duceOa the fluctuations responsible for the tail of the Et and Ezdc distributions (Eqs. (3) and 
(4) give only the average values at each b) . They have been introducedi'i in the model by multi- 
plying Ny'^^{b) in Eq. (||) by F{b) = Et / ET{b), where Et is the measured value of the transverse 
energy. Likewise, for the suppression versus Ezdc we multiply Ny"^^{b) by ha^Ezdc) /n-Aib) 
where nA{EzDc) is obtained from Eq. (4): UAiEzDc) = {AE-in - EzDc)/Ein{l - a). 

The results are presented in Fig. 2 and compared with the NA50 dataB. We see that from 
peripheral collisions up to the knee of the Et distribution, the data are well described. However, 
beyond the inflexion point at the knee, the decrease in the data is sharper than in the model. 
Note, however, that the data beyond the knee are obtained with the so-called minimum bias 
(MB) analysis. Only the ratio J /ip over MB is measured and it is multiplied by a theoretical 
ratio DY/MB. In the model used by NA50, this ratio is essentially flat beyond the knee - due to 
the fact that the tail of the Et distribution of hard (DY) and soft (MB) processes is assumed to 
be the same. The behaviour of the MB NA50 data beyond the knee can be explained 113 by the 
combined effects of a small decrease of the hadronic Et in the J/ ip event sample (due to the Et 
taken by the J/^p trigger), together with the sharp decrease of the Et distributions in this tail 
region. This phenomenon does not affect the (true) ratio J/V' over DY (obtained in the NASD 
standard analysis), but does affect the one obtained by the so-called minimum bias analysis t3. 

The resultsEJ of the model as a function of Ezdc are presented in Fig. 2 (b). The centrality 
dependence of the J/tp suppression is reasonably well reproduced except for the broad bump 
centered at Ezdc ~ 12 TeV. No such structure is seen in the Et distribution at the same b 
[Et ~ 85 GeV) . Note that the suppression beyond the knee {Ez dc ^ 5 TeV) is well reproduced 
by the model. This was to be expected since the tail of the Ezdc distribution is not affected by 
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Figure 2: a): Ratio J/^ over DY versus Et compared to NA50 datau. The full curve is the theoretical predictiem 
with Et fluctuations and the dashed line contains Et fluctuations and the Et loss induced bvithe J/ip trigge£2l. 
b): Ratip-iJ/?/' over DY versus Ezdc compared to preliminary results presented by NA5O0. The full line is 
obtainedLJ computing for each h, the value of Ezdc and the value of the ratio _R, taking into account fluctuations 
and changing the normalization by a factor 0.94. In both flgures the dotted line is the NA50 absorption model, 
fltting pA and SU. The used parameters are aabs= 4.5 mb, (Tco= 1 mb, q= 0.62, a= 0.076 and A^/= 1.15 fm~^. 



the presence of the J/ijj trigger. This provides a consistency check of our interpretation of the 
suppression versus Et beyond the knee. Note also that the data of Fig. 2 (b) for Ezdc < 28 TeV 
exhibit a J /ip suppression significantly steeper than the one in the NA50 absorption model - 
indicating that the anomalous suppression is already present in very peripheral collisions. This 
confirms the trend already observed in the suppression versus E^- 
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